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Abstract 

Symmetric hyperbolic systems of equations are explicitly constructed for a gen- 
eral class of tensor fields by considering their structure as r-fold forms. The hy- 
perbolizations depend on 2r — 1 arbitrary timelike vectors. The importance of the 
so-called "superenergy" tensors, which provide the necessary symmetric positive 
matrices, is emphasized and made explicit. Thereby, a unified treatment of many 
physical systems is achieved, as well as of the sometimes called "higher order" sys- 
tems. The characteristics of these symmetric hyperbolic systems are always physi- 
cal, and directly related to the null directions of the superenergy tensor, which are 
in particular principal null directions of the tensor field solutions. Generic energy 
estimates and inequalities are presented too. 

PACS: 04.20.Ex, 04.50.+h, 02.40.-k, 02.30.Jr 

1 Introduction 



First order symmetric hyperbolic systems of partial differential equations are of paramount 
importance in mathematical physics. Since the pioneering work by Friedrichs many 
such systems have been studied in different areas. In gravitational physics, they are 
relevant both in studying particular physical systems on given spacetime backgrounds, 
and in considering the gravitational field equations themselves — be them Einstein's field 
equations or more general possibilities — . Their importance resides on the availability of 
very powerful theorems of existence and uniqueness of solutions, under mild continuity or 
differentiability assumptions. 

Hyperbolicity of a system of partial differential equations is a concept which generally 
implies the well-posedness of the appropriate Cauchy problem, that is to say, a well-defined 
initial value formulation. The general reason behind this is the existence of norms — on 
the space of solutions — which are well-behaved under the evolution denned by the system. 
One of the purposes of this short paper is to bring out the relevance that the so-called 
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'superenergy' tensors have in the definition of these norms and on the hyperbolizations 
of general systems in Lorentzian backgrounds of arbitrary dimension. 

Symmetric hyperbolic systems for most of the known physical fields can be set up, 
in both flat and non-flat spacetimes, with more or less difficulty For an account of 
these, one can consult fT7\ [321 12] and references therein. In these cases, the background 
spacetime is a given Lorentzian manifold and the field equations can be written with 
the metric, connection and curvature variables as known data. A much more difficult 
problem is that of solving the field equations for a physical field coupled to gravity, for 
the background spacetime (metric, connection, etc.) must be built at the same time by 
solving the corresponding field equations. In this case, the causal structure inherent in the 
symmetric hyperbolic systems and that corresponding to the spacetime thus constructed 
should agree, leading to causal propagation of gravity and physical fields. The importance 
of superenergy tensors in the study of causal propagation was already made explicit in 

0HU]. 

In recent years there has been an increase of attention concerning symmetric hyper- 
bolic systems for Einstein's field equations due to the exigencies of the developing area 
known as Numerical Relativity, see e.g. jHJ |2U E3 ES] and references therein, where initial 
data are evolved by numerical integration of the appropriate field equations. As a mat- 
ter of fact, the history of symmetric hyperbolic systems in General Relativity is largely 
related to the different developments concerning the Cauchy problem for the Einstein 
field equations, see (TSJ EH] [21] and references therein. For very good recent surveys of 
the gravitational Cauchy problem, see [23 121] • A new input to this problem came with 
the work by Friedrich [221 EH] where the Bianchi identities for the Weyl curvature tensor 
were included in the system, and the Bel-Robinson 'superenergy' tensor was utilized to 
estimate its strength. This was later improved in ^B] , where only physical characteristics 
were obtained and matter sources could be included. These results relied on an 'electric- 
magnetic' decomposition of the Riemann tensor, and on its derived "superenergy", both 
introduced by Bel [3] in 4 dimensions many years ago. Improved results of this approach 
were found in pQ, where an integral inequality for the Bel superenergy was introduced 
and causal propagation of the Riemann tensor in vacuum was also obtained. In these and 
later papers [T71 ITS] , all the results were derived in four dimensions. 

An interesting approach to these results was given in 0, by just presenting direct 
hyperbolizations — in the sense of |2Z] — of the Bianchi equations. Again the results were 
obtained in 4 dimensions. However, the whole argument in [Hj is valid, mutatis mutandis, 
in arbitrary dimensions except for the crucial final step where the positivity properties of 
the system's symmetric matrix was proved. To that end, in [U] spinors were used 1 , which 
restricts the result to 4 dimensions exclusively. 

This is a key point in the present paper: for some time now it has been known that the 
electric-magnetic decomposition of the Riemann tensor, and actually of arbitrary tensors, 
as well as the superenergy construction and the positivity properties of the superenergy 
tensors, hold in arbitrary dimensional Lorentzian manifolds and furthermore, universally, 

Actually, by using spinors the hyperbolizations of the Bianchi identities, and of quite general spinor 
fields, are rather obvious, see [23 EID E] 
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that is to say, for arbitrary tensor fields. The basic references for this are jHU E2] and 
references therein. 

By using those facts, I am going to show how to construct, in general Lorentzian 
manifolds of arbitrary dimension, first-order symmetric hyperbolic systems of equations 
for general tensor fields subject to suitable field equations, and the crucial role that the 
positivity properties of the superenergy tensors play. Let (V, g) be an n-dimensional 
Lorentzian manifold with metric tensor g and signature (— , +,..., +). The archetypical 
system to be hyperbolized is 

where the unknown A w ... Mm is an arbitrary rank-m tensor field, J^...^ m = J[^ ...^ s+1 ...^ m 
and jw.-nm = J[M2...MsW+i---Mm are gi ven an d may depend on the background Lorentzian 
manifold and on the tensor field -Atm (but not on its first derivatives), and s is any 
natural number such that < s < min{m,n}. This is not the only type of system 
susceptible of study by superenergy techniques, but it will serve to make the main points. 
Other more general systems can also be treated in the same manner, and some comments 
and examples will be made in section [7| In particular, one can deal with several, possibly 
interacting, fields by either (i) letting A w ...^ m to be an inhomogeneous 'multi-tensor' as in 
|3*T] or (ii) by adding systems similar to for the other fields. The latter possibility will 
be adopted here, see subsection 17.31 The Bianchi equations with and without sources, in 
arbitrary dimension, will be dealt with in subsection 17.71 

The approach presented herein has the virtues of (i) unifying many different systems of 
equations, (ii) opening new possible applications for yet unexplored cases, (iii) providing 
a simple interpretation of the characteristics of the systems, (iv) providing in a direct 
general manner integral inequalities for the "energy" of the systems (which is nothing 
but the "superenegry density" [US]), (v) allowing us to deal with systems with only one 
of the two expressions in (JTJ) by adding 'gauge' equations, and (vi) giving, under some 
circumstances, natural divergence-free tensor fields as well as conserved quantities. 



2 The 'superenergy' construction 

First of all, some fundamental properties of the superenergy construction and the superen- 
ergy tensors need to be recalled. Everything is based on the following basic result |HH EH] , 
see also j20]: any tensor A M1 ^ m can be considered, in a precise and unique way, as an 
r-fold form, that is to say, as a tensor belonging to A™ 1 <g> . . . <g> A nr where A p is the set of 
p-forms. The number r is a well defined natural number called the form- structure number 
of ... Mm and, obviously, r < m. Similarly, the set of r natural numbers ni, . . . , n r is 
uniquely defined, and each rix is called the T-th block rank. It is trivial that Ysr=i n T = m , 
and Hx — min{n, m} for all T = 1, . . . , r hence m/n < r < m. Tensors seen in this way 
are called r-fold (ni, n r ) -forms and denoted by A[ ni ]...[ nr ,] |3*3] . A tilde on a tensor 
(^ti.../i m > or equivalently, A[ ni ]„.[ flr ]) indicates that the indices have been permutted so 
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that the first n\ indices of A M1 .. >Mm are those precisely in [ni], the next n 2 indices are those 
in [712], and so on. 2 

Once the indices are split into skew-symmetric blocks, the canonical 'Electric-Magnetic' 
(E-H in short) decomposition of the tensor, associated with any unit timelike vector u, is 
immediately obtained. Probably the most natural way to do this is by using the Hodge 
dual operators * T acting on each block PHES1I2DI- Thus, there are exactly 2 r electric- 
magnetic parts of a given tensor whose form structure number is r. 3 For each block, the 
contraction with u produces the electric part in that block, and the contraction with the 
dualized block the corresponding magnetic part. Another way of seeing this is by con- 
tracting each block with u (electric) and then taking the wedge or exterior product of that 
block with u (magnetic). Yet another even simpler way of thinking of this is by taking, in 
any orthonormal basis {e*o, e\, . . . , e„_i} with u = e*o, the Q12 ■ ■ ■ i nr -components (electric) 
and the i\ . . . ^-components (magnetic) in the selected block, where Latin small indices 
take the values 1, . . . , n — 1. All the E-H parts are spatial tensors (orthogonal to u) and 
they completely characterize the tensor field. 

Given any tensor A Ml Mm with form structure number r, its basic superenergy tensor 
i^i...A r /v{^4} has 2r indices, and it is symmetric on each of the r pairs |35j. The standard 
definition of the basic superenergy tensor is given with Hodge duals, but a very useful 
property is that they can be written, actually, independently of the dimension n and 
without any duals, see section 3 in [33]. This implies that one can actually define a 
'superenergy operator', valid in arbitrary dimension n, whose outcome when applied to 
the tensor product A ® A is the basic superenergy tensor of A. This operator reads 



ip ai...a ni ...T 1 ...r nr — 



r<j 2 ...cr. 



J XllJ.l...Xrlirpi...p ni ...vi—Vn T ~ _ -Q| U P2-Pn ± y ;/ . 

71« (^ll^ - ^l\9K») (2) 



[n r 



where 



is the Kronecker symbol of order p. In other words, the operator consists of a factor 

1 x a2—^n r ( _ _J_ STtTi 

(nr — 1)! " P " T \ ( x r9fir)pi n ^ Pl 9XrfJ.r 



2 Some simple examples are: any p-form S^...^ = S^j...^ ] is trivially a single (that is, 1-fold) p-form, 
while VyE^...^ is a double (l,p)-form. The Riemann tensor R a j3\^ is a double symmetric (2,2)-form 
and the Ricci tensor Rp^ is a double symmetric (l,l)-form. In general, any completely symmetric rank-m 
tensor is an rn-fold (1,1,. . . ,l)-form. A rank-3 tensor A a p 7 with the property A a py = —A^pa is a double 
(2,l)-fbrm and the corresponding A is clearly given by A aj g 7 = A^ a p^ = A ai p. 

3 Of course, if the original tensor has extra symmetries, or if some of its traces vanish, then there 
may be some relations between different E-H parts. Thus, for instance, the two mixed E-H parts of the 
Riemann tensor are related, and there are some extra properties for the Weyl tensor C a f3\^. The fact 
that in n — 4 there are only one electric and one magnetic part of C a p\^ is a purely dimensional result; 
for further results and explanations, see |85l 
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for each antisymmetric block [ny] of the tensor A[ ni i .[^1. The basic superenergy tensor 
of is then simply 

T^ \ \A\ = -E\ \ ff i- ff »i- T i- I "» r A hr --«,---«,, (o\ 

± \\p\...\ r p r l^M g *if'i---ArHrpi...p ni ...VL...v nr ■ rL cri...a„ 1 ...Ti...T nr - rl • K'-'J 

Observe, as remarked above, that this is independent of the dimension n. Note also that 
once this formula is known, one may let the operator act on two different tensors A^^^ 
and 5^,,,^ as long as they have the same form structure number and the same block 
ranks, thus defining an operator acting on the tensor product A® B: 

T 1 -f A Rl = — F cr 1 ...a ni ...Ti...T nr 7 r>p 1 ...p ni ...v 1 ...v nr 

A \ip 1 ...\ r p r \- n -i 1J J n I1 '\ipA--A T p,rpi...p nx ...i'x—Vn T ^cri...<7 ni ...Tl...T nr 

An equivalent way of defining this operation is (suppressing indices) 

T{A, B} = 2T {^] ~ \t{A) - l -T{B) 
so that T{A, B} = T{B, A} and T{A, A} = T{A}. 



3 Hyperbolizations 

There are three inequivalent possibilities for the system ([!} according to whether s is 
lower, equal, or greater than the first block rank n\ of A. (Of course, one could also study 
the system (JI} by antisymmetrizing over any chosen s indices of the tensor field, and not 
necessarily those including the first one. The argument will work just the same, so that 
this must be understood as a way of making the reasoning clearer without affecting the 
generality.) If s < n\ (including the case s = 0), then in fact the system is equivalent to 

^ po^n\...p m = (s + ^)Jp [p, 1 ...p ni \p ni+1 ...p m — (— 1) 1 s J[p 1 ...p 7ll ]p p, ni+1 ...p m , (4) 
so that the divergence equation V p A pM2 Mm = jp 2 ...p m follows from (JU) and thus 

jp,2—Pm ~ ( S ^)J P [PfJ-2---Hn 1 ]Hn 1 + l---Llm ~ ( _ ^ ) S ^[PP2 ■ ■ -Mnj ] ^Mnj +1 • • -Pm 

for consistency. Equations (jJJ) are very simple and easily hyperbolized — see footnote 
El below — . If s > n\ then there are equations only for the new tensor field A' Mm = 
^[p 1 ...p s ]p s+1 ...p m which now has s = n[. Finally, the most interesting case is precisely when 
s — iii, which includes the previous one. 4 

In order to find the hyperbolizations of (|1J) in this case, recall a very important property 
of the basic superenergy tensors: the tensor Q is the essentially unique (i.e., unique up 
to index permutations; see section 5 of [SHj) tensor quadratic on A m ...^ m which satisfies 

4 The system Q for s — n% can be written in an intrinsic manner by using the exterior differential 
and co-differential acting on the blocks of ?'-fold forms, see (20) and references therein. With the notation 
introduced in (20j, the system reads simply dn^A — (s + 1)J and SmA = —j. 
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the dominant property [7] . A convenient way of stating this fundamental property for 
the purposes of this paper is that, for any set of arbitrary timelike future-directed vectors 
{4\ <\..., uf,v?}, 

TW.^^KX 1 • • • u)rv^ > (5) 
(for non-zero A M1 .. -Atm ). An equivalent statement is that the vectors 

i 5a = -n,..v r {^K 1 -^Y 

are causal and future directed. It follows that the vector- valued matrices (endomorphisms 
acting on the set of r-fold (n 1; . . . , n r )-forms) defined, for arbitrary timelike future-directed 
vectors {v^ 1 , U2 2 , vtj: 2 , . . . , u^ r , vp r }, by 

/10ffl...(T„ Tl...T„ r _ M [a 1 ...a ni ]...[T 1 ...T„ r ] Ml \ 2 u 2 ?/ V ? ,Mr ( f.\ 

V pi ...p ni ...vi...u„ r — 111 t^l^2^2-XrtJ.r[p 1 ...p ni }...[u 1 ...u nr } U U 2 U 2 ■ ■ ■ U r U r \P) 

are appropriate candidates for the symmetric positive-definite matrix of a symmetric hy- 
perbolic form of (JTJ). Indeed, for arbitrary r-fold (ni, . . . , n r )-forms A[ ni ]...[ nr ] and -B[ ni ]...[n r ] 
one has that 

D a ( A R\ = D a CT 1-" CT "1 — Tl—Tnr 4 f>pl»-Pnx—Vl—Vnr — 

V — V Pl ...p ni ...V 1 ...V nr ^Vl---<Jn 1 ---Tl...T nr l-> 

= T% lA2M2 ... ArMr {A, B}v^u?v? . . . uJ-<- 

so that on the one hand 

Q a (A,B)=Q a (B,A), 

and on the other hand, for any timelike future-directed 1-form u a , u a Q a (-, •) is positive 
definite as follows from (v M = v±) 

u a Q a (A, A) = -u a P a > 0. 

Now it is easy to find hyperbolizations (in the sense of j2Zj) of the system (JTJ) for s = n\. 
First, rewrite (JTJ as 

ii/ °uo...u, > /_ , 1 M "unui ...u„ _ i°/3i.../3„„ i l, w...i/». v oA; 



y ( s _X)| y tt2-Ms ' ^ _|_ "nom-Hs J n2 \Pl-Pn 2 n ^ u u 1 ...u nr w a^a 1 ...a s p 1 ...p n2 ...T 1 ...T nr 

= (j,,. U-> no • • • ) (7) 

and now contract this with 

1 ' 1 % ^:^-i^i,[wj^]]x 



2(s-l)\ V (s-1)! 171 72 - 7 - s\s 



71 •••7s 



-2 x 



(n 2 - 1)! V 2 W2[ei ^..-n 2 ] +^2 «2[ £1 <) e2 ..,„ 2 ] — ^...e n ^ 



The result is 



V 71. ..7s ei...e n2 ...Ci— C«r a-^a 1 ...a s p 1 ...p n2 ...Ti...r„ r — c/7 1 ...7 s ei...en 2 -Cl-Cn r 

where J 71 ... 7s ei...e„ 2 ...Ci-Cn r i s the contraction of (jHI) with the righthand side of (J7J). This 
form is manifestly symmetric hyperbolic. Taking into account that the part of Q which 
depends on u 2 , v 2 , ■ ■ ■ , u r , v r is non-degenerate when acting on (r — l)-fold (n 2 , . . . , n r )- 
forms, this system can in fact be written in the simple form 

(s + l)^V [p i 7l ... 7s]ei ... en2 ... Cl ... Cnr + (-l) s sw [7l V p i 72 ... 7a]pei ... £n2 ... Cl ... Cnr = 

(a + l)^Jp 7l ... 7sei ... en2 ... Cl ...c„ r + s% 1 j 7 2...7^i...en 2 ...fi...Cn r • ( 9 ) 

Let us finally remark that the hyperbolicity of the general system (JTJ is related to the 
existence of a wave equation for v4 Ml ... Pm . This is immediate from the definition [20] of the 
de Rham operator A(i) acting on the first block, which is given — using the notation of 
footnote |U — by Am = ^(i)^(i) + ^(i)^(i)- Thus, from (JTJ) one deduces 

A {1) A = -d m j + (s + l)5 {1) J 

which reads, in index notation, 



X7PX7 A — q/? r AP i -l- S ^ S ^ /? r A~P a 1 

V v^j,,,^ Olhp^rx f_ l2 ...f_ ls \fi s+1 ...p m T ^ 1L pv[PlP2^- /Lt3— MsjMs+1— Mm 

m 

+ S RF ^u^A n2-»Ha]lis+l—fH-lpfH+l—Hm ~ S ^lPl jn2---fls]/J.s+l—fim + ( S + 1) V P Jp pi ... pm . 

1=5+1 



4 Constraint equations and integrability conditions 

In order to prove the equivalence of the hyperbolized system Q with the original ([TJ, let 
us do some counting: the number of unknowns is MC nyS where Af stands for C„ in2 • • • C n . nr 
and C n>s = n\/{n — s)\s\. The number of equations in (JTJ is jV(C„ )S _i + C„ )S+ i). There 
are more equations than unknowns, so that the system is overdetermined. However there 
are constraint equations — see [27] — which are easily computed. These are given by, for 
any (n — l)-dimensional hypersurface £ with normal 1-form iV M , 

N P -f4a-p 7 3... 7s ei...e n , 2 ...Ci—Cn r — j fr/3.---y s ei--^n 2 — Cl—Cnr ) ~~ 0, (10) 

A^[ CT V 70 A 7l ... 7s ] ei ... £n2 ... Cl ... Cnr - iV[ (J J 7071 ... 7s ] ei ... en2 ...^ 1 ...c nr . = . (11) 

Notice that the first of these is absent if s — n\ < 1, and similarly for the second if s = 
ni > n— 1. The total number of constraint equations is therefore A/"(C n _i )S _ 2 + C n _i )S+2 ) 
and given that 

C*n-l,s-2 + Cn-l,s+2 = C*n,s-1 + Cn,s+1 — C*n,s 
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the constraints are complete. 5 This leads to the equivalence of the new system Q to 
the original one (JTJ) together with the constraint equations ()10lll|) referred to a spacelike 
hypersurface with timelike normal N p . Observe that the system is causal in the sense that 
every 1-form u a such that u a Q a (-, •) is positive definite is timelike. Hence, one can define 
a well-posed Cauchy problem (or initial value formulation) of the system (JTJ) by giving 
initial data, subject to satisfying the constraints, on any spacelike hypersurface. 

The integrability of the constraint equations is not ensured, and depends on the 
Lorentzian manifold and on the righthand sides of The integrability conditions are 
ruled by the squares of the operators , 6^ — see footnote 0] — , which are in turn 
governed by the Ricci identity. Using formulas (23) and (24) in (201 these integrability 
conditions are 

Z i=s+l 
m 

pp3...p m ~ n ^ Mi-V 1 ^ P3...p 3 Ps + l---Pi-lPPi + l---Pm • 

Z i=s+l 

Therefore, in flat spacetimes these are rather simple (and easily satisfied), but they may 
set very strong restrictions on the form of J and j depending on the algebraic properties of 
the Riemann tensor. The exceptional cases are (i) when ni < 1, in which case the second 
condition is absent; (ii) the case n\ > n — 1, so that the first condition is trivial; and (iii) 
when the unknown A has r = 1, i.e., it is a single m-form and thus n\ = m. Then the 
previous integrability conditions are simply V[a Jpp 1 ...p m ] = an d ^ p jpp 3 ...p m = 0, which 
can be enforced in any Lorentzian manifold — for instance by setting J = dK and j = 5k 
for arbitrary m-forms K, k — ; see subsection 17.11 

These integrability conditions are reminiscent of the well-known Buchdahl compati- 
bility conditions which arise for spin greater than one in four dimensions ^21 El , see also 
the discussions in jSUl El ■ Of course, they may drastically reduce the number of solutions, 
or even forbid their existence completely. As a matter of fact, the system contains the 
natural generalization of spin-S* fields for any value of S to arbitrary dimension n. 



5 Characteristics 

Let us now consider the important question of the characteristics of the system ([!}. By 
definition the characteristics are given by the directions c a such that there are non-trivial 
solutions A M[n2 ]...[ nr ] of 

f\OL a 1 ...a a Pl---Pn2--- T l--- T ^r I _ r> 

L QV 71. .. 7s e 1 ...e„ 2 ...C 1 ...Cu r ^l-<J s pi-Pn 2 -ri...T„ r U. 

5 By the way, the very same Q a pl'.'.'. P 2l'/.'.ul'.'.'.u2l as defined in © hyperbolizes directly the system 
The resulting system is in fact exactly @. However, in this case the constraint equations are not i|l(Jlll[l 
but simply A r [aV 7o ]4 71 ... 7 , ei ...e„ 2 ...Ci...C„ r = A r [ o -/ 70 ]7i...7»ei...en 2 -Ci-C» r and tne y are a S ain complete. 
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This condition can be proven to be equivalent to (setting {Q} — e\ . . . e n2 Ci ■ • • Cn r 

to alleviate the notation) 

sw p c [7l i p 72 ... 7s]{ n} + sc p v hl A p ~ /2 ... ls]{n} - v p c p A 11 ,„ ls{n} = (12) 

whose more general solution is given by either 

i. the null 1-forms c p such that c p Apj 2 _ Js ^u} = and C[ 7o A n ... 7s ]{n} = 0, so that the 
solutions have the form A lx ... 7a {m = c [-aX l2 ... ls ]{n} with c p x p ... 7s {n} = 0. 

ii. the spacelike 1-forms c p such that v p c p = so that the solutions take the form 
^7i...7 S {Q} = W [7i c 7 2 l/73...7 S ]{^} where, without loss of generality, y l3 ...~, 3 {n} can be taken 
orthogonal to both c and v in the first block. 

As a first conclusion, the characteristics are all physical, in the sense that they define null, 
or timelike, hypersurfaces of propagation. A more interesting conclusion is that these 
characteristics can be characterized, in all cases, by well-defined null directions associated 
with the solutions. In order to prove this, let us define for the case EH two null directions 
{£+,£'} such that Span{c,w} = Span{/+,/-} and £ +p £- ti = -1. Then 

W^V.^IM = o, i;£^A% 2 ^ s]m = 0, 

which implies that both £~} are in particular principal null directions (see e.g. Defini- 
tion 2 in [SI]) of the tensor A s ir n2 i r n j. The same can be said of the null c in the first case 
above. It is known (31] that k is a principal null direction of A s ]r n2 i r^i if and only if it 
defines a principal direction — see Definition A. 2 in [SB] — of its basic superenergy tensor 
that is to say, if and only if k Xl k pl ■ ■ ■ k Xr k^T Xl ^...x r ^{A} = 0. Due to 
the general properties of causal tensors 01213 ; these principal directions are necessarily 
null. Furthermore, in the case under consideration one is dealing with a special type of 
principal directions, because they are related to the first block exclusively. This implies 
that the characteristics of the system Q are defined by the — necessarily null — k p such 
that 

k^k^T XlPl ... Xr , r {A} = 

where k p stands for either c p in case or £ ±p in case EH 

In summary, all characteristics of the system (QJ) are physical and they are directly 
related to special principal null directions of the corresponding solutions. 

6 General energy estimates and inequalities 

The hyperbolizations defined by (jUJ) depend on 2r — 1 arbitrary timelike future-directed 
vectors (thus, on one timelike vector for m-forms, on three timelike vectors for double 
forms such as the Riemann tensor, and so on). Selecting all of these vectors identical and 
thus equal to a unique (unit by convenience) timelike vector, say u, then u a Q a (A, A) is 
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exactly the "superenergy" density 6 of the tensor A Mi ... Mto . This is denoted by Wa(u), it 
can be defined alternatively as 

W A (u) = T Xl ^... Xr ^{A}u Xl u^ . . . u x ^ 

or as 

z \ x=1 n r ij 

where 

h^ u (u) = g^v + 2u IM u u , 

and it is equal to half the sum of the positive squares of all the electric-magnetic parts 
of the tensor A M1 .. >/lm [55] . The supernergy density is also equal to half the sum of the 
squares of all the independent components of A M1 ... Mm in any orthonormal basis {e^} with 
e*o = m (this is the typical mathematical 'energy' of the tensor A fll ___ IMn ): 

i / r i \ n-1 

^A(eo) = T 00 ... {A} = - JI -7 E l^,..,J 2 - 

This function allows us to obtain estimates because the dominant property implies [35] 
that Wa{u) = T o...o{^4} > |^)ti...^ m | in any given orthonormal basis. One can also obtain 
integral (in)equalities related to the system (JTJ). As a matter of fact, this quantitiy and the 
properties of the full superenergy tensors were already used to prove the causal propaga- 
tion of the gravitational field in (TJIj, and of general physical fields in arbitrary dimension in 
. (The causal propagation along characterisitics can also be deduced from the study of 
the propagation of discontinuities; see [23 E3] section 7.3 in [35], and references therein.) 



Of course, the Bel-Robinson and Bel tensors in four dimensions |U [5] E] have been re- 
peatedly used to obtain estimates [2E] > and to prove the hyperbolicity of Einstein-Bianchi 
equations [221 1231 HI EH El UHj and the non-linear stability of flat spacetime |XSJ . 
Generic integral inequalities are derived by first defining EI] 

w(t) = / T Xlfil ... XrlM .{A}N^N^ . . . N K N^ 77 = 

JD+(E t0 )nJ-(s t ) 

J" ^ T Xllll ... Xrllr {A}N^ . ..N^N*~da x ^ t ^j dt' > 

where £f is a compact achronal set (usually a — piece of — a hypersurface) , D + (E to ) 
its future domain of dependence and r\ denotes the volume element n-form. Due to its 

6 This name may be inappropriate sometimes, because this is the traditional energy in relevant cases 
such as the electromagnetic 2-form, or in the case of the massless scalar field <j> by using dip as the tensor A. 
However, is some other cases the relevant quantity has been traditionally called "superenergy" , starting 
with the Bel- Robinson definition for the Weyl tensor 0] or the Bel one for the Riemann tensor ISIE], and 
also especially for higher order superenergies such as the one defined by using the covariant derivative of 
the electromagnetic 2-form, see ^JES], or the second derivative of the scalar field, see |35], references 
therein, and subsection 17. 71 
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global hyperbolicity, D + (T, to ) is foliated by spacelike hypersurfaces T, t = {t = constant} 
with future normal N = — dt whose volume element (n — l)-form pointing along N is 
denoted by da x \Y, t . Observe that w(t) = only if T\ lfll „.\ rfJir {A} = 0, which is equivalent 
to A^...^ = 0, on D + {T lto ) PI J~(£ t ). It is very simple to see that dw(t)/dt > and 
moreover, using the Gauss theorem one derives for some constant M jH] 

< Mw(t) + jf T Xl ^... Xr ^,{A}N^ . . . N^N^da* | StQ 

+ / V p T p ^.. Xrt , r {A}N^ . . . N Xr N^ r r] (13) 

ij-(Et)nZ?+(S t0 ) 

so that the divergence of the superenergy tensor controls the growth of w(t). But this 
divergence can be obtained easily using (J3J and (J2J) 

X7 TP , , / 41 — , (J 1 ...(J ni ...T 1 ...T nr 

V p-t fll\2^2---^r^r\^ L I Al W ' ' ' A ' ' pr Pi • • -P«l • • ■ - M ™r 

Jf)2---Pn 1 1 ...a n2 ...Tl...T„ r ^p,l 

~J ^ L ^lP2---Pn 1 cri...CT n2 ...T\...T nr 

Ul _ f 7 A/}\...p ni pi...pn 2 ..Mi...v nr 

\ U /J,lf3l...l3 ni <Tl...crn 2 ...T\...Tn r - n - 

I 7 Pl.../3n 1 pl...pn 2 ...vu..v nr A \ 
~ u p,l snpi...fj ni cri...o-„ 2 ...Ti...T„ r J ■ 

Notice that a sufficient condition for the superenergy tensor to be divergence-free is that 
j = and J = 0, see also [31]. In this particular case, for instance, as well as in the 
variety of cases in which the last integral in (|T5j) vanishes, (jTBjl can be easily resolved by 
using the Gronwall Lemma so that 



w 



(«) < ± (e Mt - l) J T Xl ^ Xrflr {A}N^ . . . N^N^da^ | KtQ . 



Similar arguments can be used whenever j and J are polinomic on the tensor A Ml ... Atm , or 
even in more general cases. 

Observe finally that Tx lfll ...\ rflr {A}N 111 . . .N^N^ is conserved — that is to say, it is 
divergence-free — if V P T P f j, 1 \ 2P2 ... XrfJ , r {A} = and N is a Killing vector [33] . 

7 Examples and applications 

In this section some selected applications and relevant cases of the general system (JTJ) are 
presented. 

7.1 m-forms 

The simplest application of the results is to the case with r = 1. The equations are simply 

dA = (m + 1) J, 5A = -j 



11 



which are Maxwell-like equations for ordinary m-forms A^...^ = A^...^]. These are of 
interest in higher dimensional theories such as string theory, supergravity, etcetera. Of 
course, they include the electromagnetic field equations in any curved background of any 
dimension by setting m = 2, J = and letting j to be independent of the 2- form A. 
From the previous equations one also gets a wave equation 

AA = -dj + (m + 1)5 J 

where A is the de Rham operator. 

As remarked at the end of section 0] the integrability conditions are satisfied whenever 
dJ = and 5j = 0. The symmetric hyperbolic systems Q read in this case 

(m + l)w p V [p A 7l ... 7m] + (-l) m mt» [71 V p A 72 ... 7m]p = (m + l)v p J ni ... lm + mt» [71 j 72 ... 7m] 

for any timelike vector field v, which in fact provide all possible hyperbolizations. They 
are determined by 

Tx,{A,B} = —L^ (A { ^B, )p2 ... Pm - ^ PlP2 --- pm B PlP2 ... Pm ) 

which defines the standard energy-momentum tensor of A by Tx^A} = Tx^{A, A}. 

Observe that the characteristics of this system are determined by the null eigenvectors 
ofT AM {A} 0Eq. 

7.2 Extra algebraic relations 

For the general case of r-fold forms with r > 1, apart from the antisymmetries of each 
of the r blocks one may have to deal with some additional symmetries involving indices 
of different blocks. A typical example is that of the Riemann tensor, which is a double 
(2,2)-form satisfying the first Bianchi identity R a [pxfi] — which implies the symmetry 
between the two pairs of indices. Trace-free properties are also in this category. Another 
typical possibility, in the case of having several fields, is an algebraic relation between 
them. An example of this is the second Bianchi identity with sources for the gravitational 
field; see below in subsection 17.71 

The optimal way to deal with these algebraic properties is to simply ignore them on 
the unknowns, so that the tensor A is not assumed to satisfy them, and add them as initial 
conditions of the Cauchy problem defined by the system of equations. If these properties 
are preserved by the evolution, then the corresponding solutions will satisfy them. 

The strategy, schematically, is as follows. Suppose that A is an r-fold (ni, n 2 , . . . , n r )- 
form, and assume that there is an extra symmetry property mixing (say) the first and 
second block, such as 

A[ fJll ...fi niV1 ...Vg\u q+ \...Vn 2 Tl...T nr ~ 0- 

Define the r-fold (ni + q, n% — q, . . . , n r )-form 

P = A, i 

1 [ll... P „ 1 + q Vl...V„ 2 - q ...Tl...T nr ^■[fll... P „ lPni +l...fl ni +q\Ul...U„ 2 - q ...Tl...T nr 
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and derive differential equations for this tensor field from the original system (0). For 
example, 



or 



V ' [fi Pf J , 1 ...fi ni+q ]i/ 1 ...u n2 ^ q ...Ti...T rLr J[fM)fMl--.tln 1 thi 1 +l»'IM 1 +q]vi..Mn 2 -q.--Ti...Tn r 



X7 p P - (1 - n 

V 1 pH2---L-t-n 1 +qVl..M n2 - q ...Tl...Tn r ^ Til + CJ- l,q J f!2 . . ./inj +q V\ . . M n . 



2 — q...T\...Tn r 

+ ( — 1) 1 q {n\ + l)C ni -|- g _l ig _l J P [fJ, 1 ... f J, ni+q ]pV2...V n2 - q ...Tl...T nr 
V 1 ) ">l^ni+q-l,q-l v [pi^ p2---Pn 1 +q]p^2---l'n2-q---ri---rn r ■ 

The goal would be to get an homogeneous system whose unique solution, given initial 
conditions such as the vanishing of P, is the zero solution. In this particular example a 
necessary condition is, for instance, J[ Mm ...^ ni+1 ... Mni+9 ]^...^ 2 _ g ... Tl ... Tnr = 0. 



7.3 Mixed, or interacting, systems 

As explained in the Introduction, one can deal with as many fields satisfying as desired, 
and the righthand sides may then depend on all of them. For example, consider the system 
(JH) together with 

^ Ipo-^fii ...fi q ]p q +i.. .p p — Jpo---ppi ^ P A[p(j, 2 ...fi q ]fi q+ i ...p p — jp2---pp (14) 

where now J flQ ... flm ,j fl2 ...^ tm , ^...^3^...^ may in fact depend on both tensor fields A^...^ 
and A^...^. Of course, this can be done for as many fields as desired. 

The hyperbolization of this mixed system is achieved by just hyperbolizing each of the 
systems separately as shown in sectional so that the relevant tensors are the corresponding 
superenergy tensors of A and A. If f is the form structure number of A, then they depend 
in general on 2(r + f — 1) arbirary timelike vectors. The total superenergy density — and 
hence the total function w(t) — is just the sum of the respective superenergy densities of 
the tensor fields A and A. Note, however, that the corresponding superenergy tensors 
have different ranks, 2r and If respectively, if r ^ f. Thus, in these cases one simply 
cannot add them tensorially. Nevertheless, when r = f this can be done and in fact this 
happens in many cases of physical relevance; see subsection 17.71 

In the case of a mixed system (|T|)- ()14j) . there may be some algebraic relations between 
the different fields, such as 

F L {A,A,x) = 0, L = l,...,l. 

These relations should be treated as extra algebraic relations in the manner explained in 
subsection 17.21 and thus assumed as initial conditions, checking their preservation by the 
evolution of the system. 
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7.4 Rank-1 blocks 



Now, let us consider the important case in which the tensor field A^,..^ has one index 
with no antisymmetry properties with any other index. In other words, there is a block, 
say the first, with n\ = 1. This always includes 1-forms and fully symmetric tensors. It 
will also be useful in the study of 'higher order' systems, see subsection 17. 51 Particularizing 
((TJ) to this case the system reads 

where J^m-nm = J\jMom]-i*m- This s y stem is always hyperbolizable as 

2v p V[pA^ 1 ]^ 2 .. ^ m — v /I1 'V f> A pll2flrn = 2v p J plllflm + v iii—iXm 
for any timelike vector v and this is equivalent to (fT5*j) if the constraint equations 

are added, where N is the normal to any (n — l)-dimensional hypersurface. These con- 
straints are integrable only if 



Z i=2 



■Un 



Observe that j^ 2 ...^ m is missing in the constraint equations and their integrability con- 
ditions, and thus it is not restricted in any way. This has important implications. For 
instance, suppose that only the first equation in (|15|) is given. One can obtain a symmetric 
hyperbolic system in this case by simply supplementing the second equation in (|T5|) for an 
arbitrary tensor field j At2 ... Atm . In other words, the divergence on the first index of .<4 Ml ... Mm 
can be considered as a gauge, not appearing in the original equations, and therefore it can 
be prescribed at will. The resulting system is always hyperbolizable as shown in section 
EH and one can find the solutions to an initial value problem either depending on j M2 ... Mm , 
or alternatively, for any given particular explicit j M2 ... Mm . 

This type of gauge equations can be considered in general, not only for rank-1 blocks; 
but then the supplementary equations are subject to the corresponding integrability con- 
ditions. This general case will be commented upon in subsection 17.61 



7.5 Higher order systems 

In this subsection two kinds of 'higher-order' systems will be considered. First, a standard 
and rather obvious construction of a symmetric hyperbolic system for the first derivative 
of any tensor field is presented. Then, the question of higher-order partial differential 
equations is analyzed using the typical procedure to rewrite them as a first-order system. 
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7.5.1 The elementary higher-order systems 

Given any tensor field -ff Atl ... Atm _ 1 , subject or not to any differential equations, a symmetric 
hyperbolic system of equations for its covariant derivative 

V^-H^a— Mm = ^Ml— Mm (16) 

can always be built. If the original H^...^^ is an r-fold (m, . . . , n r )-form, then A^...^ 
is an (r + l)-fold (1, n\, . . . , n r )-form. In particular, the rank of its first block is 1 and this 
block has no (anti)-symmetries in general with any other index. Furthermore, the Ricci 
identity implies that 

Y m 

^[MoA«l]M2— Mm = ~T. ^ R P munu-i ^M2-- IM-lPH+l ■■■f J m ' (X^) 
Z 1=2 

Therefore, the situation is that described in the previous subsection 17.41 By adding the 
gauge equations 

V p Apfj, 2 _fj, m — j^ 2 ...^ m (18) 

for arbitrarily chosen j^ 2 ...^ m — depending on the background, A and H — , one obtains a 
first-order mixed symmetric hyperbolic system constituted by the equations (jlb117H8|) . I 
will call this system the "basic higher-order system" for H . (Observe, by the way, that 
the last equation (fTHj) is equivalent to a wave equation V p V p if Ai2 ... Atm = i/^-Mm)- The 
basic higher-order system is mixed in the sense of subsection 17.31 and the corresponding 
parts are given by (i) Eq. ffTBj) by itself and (ii) the two equations (jl7H8j) for v4 Ml ... Atm . The 
part (i) is of type (jlj) and thus hyperbolized in the way outlined in footnote |5] Besides, 
its integrability conditions are always satisfied by virtue of (|17j) . The second part (ii) is 
of the generic type (JTJ) treated in this paper, and belongs to the particular case discussed 
in the previous subsection. Its integrability conditions are also satisfied by virtue of the 
second Bianchi identity for the Riemann tensor and (fTo^ itself. 

A more drastic way of building a symmetric hyperbolic system is to give the whole 
symmetric derivative of A Atl ... Atm , and not only its trace. Thus, instead of (|T%|) one can add 

V( At0 v4 Atl ) p2 ... Atm = X/LtoMl/te— Mm (19) 

for arbitrary x W)MlP2 ...^ m = X( w >mi)m2...aw Tnis leads > together with (JTTj), to a system of 
type flU). 

Of course, if l was subject to some differential equations, then the choice of 

j/i 2 ...Mm m (H3) or XMoAtiM2 - Mm i n (f*^ must be done accordingly, in a compatible manner. 
From the perspective of the basic higher-order system, any differential equations on H 
can be seen as extra algebraic relations in the sense of subsection 17.21 and they simply 
read 

F L (A,H,x) = 0, L=l,...,l 

for some functions Fl- Thus, as a first criterion for the solvability of any equations on 
H fll ...^ m _ 1 , one can always construct its basic higher-order system, hyperbolize it, and 
check if the original equations, considered now as initial conditions of the previous type, 
are preserved by the evolution of the system. 
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7.5.2 Higher order partial differential equations 

Consider now the case of a system of higher-order partial differential equations for an 
unknown tensor field L M1 ... Mm . Assume that the order of the system is k so that it can be 
written as 

/<n pl -%,...^V P] . . . V Pk L^-^ = I n (20) 

for some functions In, Kn pl '" Ph m...^ m depending on the background Lorentzian manifold, 
and on L m ...^ m and its covariant derivatives up to order k — 1. Without loss of generality 
assume that i^n ' /n---Mm = ^ci pl '" Pk tn...n m - By introducing auxiliary fields correspond- 
ing to each of the symmetrized derivatives of L fM1 _ . -Atm up to order k — 1 as follows |2Tj 

A U2 ... Ui )fn. ..fim — -Au 1 ...u i pL 1 ...iJ, m t i — 2, . . . , k — 1 

one can then construct k — 1 basic higher-order systems in the sense above by using the 
corresponding Ricci identities when needed: 

^ m 

V^A^^^j. ..^ m = Aj^ V2U3HX... /Mm + R P ^jUz{u2^Vi)^i...^ m ~ 7T -^Vl ^2>3 ^PMl ■••Mm 

+ 77 J]] V\R P lMjV3V2 "I" ^ V2-R P LljVAVlj -^l.../ij_ip/^ + i.../i m (23) 



and so on for all z = 1, . . . , k — 2, together with the last equation derivable from a Ricci 
identity, namely 

V^Aj/jjj/g... UkfMl— Mm — •f~V\WxVz...VKil\...ll m (24) 

where ^^^...^^...^ = .F^,,,,]^...^)^...^ is a known function depending on the back- 
ground spacetime and the A vl „. Viltl „. llm for z = 1, . . . , fc— 2. The result is a mixed first-order 
system formed by (j21H24)) together with the original higher-order equation fl2()D rewritten 

as 

Kn pl Pfc ^i...^ m V pi Ap 2 ... Pkfll ,„n m = In (25) 

where In, Kn pi "' Pk ^...^ are now considered as functions of £ m ... Pm and A^...^...^ for 
i = l,...,k — 1. Each of the systems (JZH), (|2*2*j) . (J2HJ), etcetera, are of the type (JIJ), 
and thus easily hyperbolized as indicated in the footnote El The constraint equations 

for these systems are complete and integrable by virtue of ()22I23[ J2I|, respectively. 

The remaining system is formed by (|2*%jl and (|25[). for which the constraints are complete 
and integrable too. The question of whether this system is hyperbolizable or not depends 
on the explicit form of Kq P1, " Pk p 1 ...p m - Some simple hyperbolizable examples leading to 
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symmetric hyperbolic systems obtained by the procedure explained in section El are (see 
also j23) 

TSPl-Pk - n Pl(p2XP3 . . . XPk) \Pl-Pm 

ZV H pi---p m if u a Li u a k ^^...fim 



or 



where in both cases A™'"^™ represents the product 5^ • ■ ■ <5^™ with the appropriate (anti)- 
symmetrizations on the indices H\ . . . \i m according to the r-fold form structure of L^...^ 
or the symmetries between its indices. 

Of course, a particular explicit very simple example of the above is the massless scalar 
field (f) which satisfies the second order equation 

VVp0 = . 

Then, one can construt a first order system on the variables (0, M ) by means of 

V M = </v, V[„0 M] = 0, V P P = 0. 

This is the case of a 1-form, m = 1 in subsection 17.11 The hyperbolizations are defined 
by the energy-momentum tensor of the scalar field 

T X/t { V0} = - \gx^ p <\> P . (26) 

Thus, the characteristics of this system are determined by the null eigenvectors of T\ p { V0} 
[3 EE!, an d they depend on the causal character of M . 

7.6 Adding 'gauge ' equations 

Consider the case when only the first part of the system (0) is given, namely 

^ \pO^Pl---Ps]Ps+l---Pm ~ JfJ.0---P.rn 

and assume, without loss of generality, that s = n\ is the first block rank of A seen as an 
r-fold form. These equations by themselves are underdetermined if s > (n — 1)/2, because 
there are NC ntS unknowns and only NC n ^ s+ i equations. Therefore, for these values of s, 
there are no possible hyperbolizations of these equations as they stand |27j . 

However, a feasible procedure for general values of s is the following. Observing that 
the symmetric part (with respect to an index in the first block) of the covariant derivative 
of A is not restricted in any way by the given equations, one can add these, fully or partly, 
arbitrarily as 'gauge' equations. Thus, for instance, one could add either of 



^ ( P 0^-pi)P2---PsPs + l---Pm Jpo-.-Pml ^ ^ ^pH2—Pm jp 



2---Pn 



for arbitrary J^^^... ^ — J^m)^— n s ]fj, 3 +i— t*m or jp2---p s ps+i---p m ~ j[p2---ps]p s +i---p m - This 
leads, in the first case, to a system of type and in the second to a system of type ((H). 
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Any of these completed systems can be hyperbolized as explained in previous sections. The 
only restrictions to be placed on j A10jU1/i2 .../Lt m or j w ... Msft+1 .. are the necessary integrability 
conditions written down in section Then, one can find the solutions to an initial 
value problem depending on the presribed j P2 ... Pm (or J popiP2 ... Pm ). Alternatively, one can 
solve the system, depending on initial data, for any given particular explicit j P2 ... Pm (or 

Similarly, if only the second in (0) is given, one can add the first in Q as a gauge for 
arbitrary J P0 ... Pm , subject to its integrability conditions. 

7.7 Bianchi identities with and without sources in arbitrary di- 
mension 

As an important physical application of all of the above, let us consider the so-called 
higher- order field equations for the gravitational field, see [21], which are essentially the 
Bianchi identities for the Riemann tensor. Letting K a p\ p = K[ a /3][A;u] be any double (2,2)- 
form, not necessarily satisfying any extra symmetry property such as K[ a px]n = 0, the 
'Bianchi' equations are 

V [7 K a/3]x „ = 0. (27) 

This is an underdetermined system when n = 4, but not for higher values of the dimension 
n. In any case, as explained in the previous subsection 17.61 one can add for arbitrary n a 
'gauge' part 

V p K pPXp = j pXll (28) 
where j pXpj = j p [ Xp ] is restricted to satisfy the integrability conditions 

V p JpXp = l - (R p Xtv K t \ p + R% TV K T \ P ) (29) 

but is otherwise arbitrary. 

These equations can be attacked in two completely different ways, both from the 
mathematical and physical points of view. The first possibility is taking the background 
Lorentzian manifold as given — and therefore its Riemann tensor is known data — , and 
treating K a p Xp as an unknown tensor field to be determined. This has some interest from 

(i) the mathematical point of view, as one of the solutions is the Riemann tensor itself and 

(ii) from the physical viewpoint, because it serves the purpose of studying "spin-2" fields 
on given spacetimes, such as a linearized field in flat spacetime and similar cases. The 
second possibility, which is much more difficult to handle and has a richer structure and 
more physical interest, is the case when the background spacetime is not given, but one 
wants to find the metric tensor field and corresponding connection alongside the solutions 
for the higher-order equations (J27|) . In this case, one has to build a hyperbolic reduction of 
the Einstein's field equations for some matter content where the unknowns contain metric 
components, connection variables, and curvature components. The way to do this in four 
dimensions, which can be translated in principle to higher dimensions, can be consulted 

in |221 ESDI DEI El DEI- 
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In both cases, the hyperbolization of the Bianchi equations by themselves, or as part 
of a larger mixed system, should be achieved. Considering the system (j27H28j) . and by 
the general results of this paper, the constraint equations are complete, a symmetric 
hyperbolization is provided by the superenergy tensor of K a/3 x p (its Bel tensor, see HH 

my- 

T a p X p {-^[2] [2]} = K ap XaKi3 P ^ + K ap ^ a Kp P \ J 

-lg a/ 3K pTXa K p \° - \g^K apaT K^ T + \g amp K pTav K pT ™ (30) 

2 Z 8 



and the characerisitics of the system (J27H28j) are determined by the principal directions 
of T Ql3Xfl {K[2][2]} such that 

i a e^T a ^ p {K m] } = 
which are the principal null directions of the corresponding i^[2][2]> defined by 

e% a Kp] P x» = o. 

Furthermore, the divergence of T a p\ p |i^[2][2]} can be easily computed usign (|27H28|) to 
give 

V a T al 3\ p {i^[2][2]} = jp\<jKp P \j° + jp^K^x 7 - -j pTa K f3 pTa . 

and this rules, as explained in section EJ the (super)-energy estimates. 

All of the above can be said in general. Nevertheless, one can consider the important 
particular case of Einstein's field equations 



— c 



or similar problems for most general equations relating the Ricci tensor R pv and the 
energy-momentum tensor T pv of the gravitational sources, which must be divergence-free 
by assumption 

VT W = 0. (31) 
Given that the Riemann tensor satisfies 

Ra[/3\p] — =>- R a /3Xp, = R\[iaf3 (32) 

from (|77j) one immediately deduces 

VRppxp = ^[\Rp]0 
so that an appropriate choice for the arbitrary jp X p, of ((23) is 

3p\n = (V[aT m]/3 - ^^9/3[p^x]T^j , T = T p p . 
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This choice is needed if one wishes that the solutions of (j27H28j) include the Riemann 
tensor of the Lorentzian metric. The integrability condition (j2Uj) reduces now, due to (p?T]) 
which implies *V p j p xp = 0, to 

R p XrvK TV pll + R% TV K™ Xp = (33) 
while the integrability conditions of the original (}2Tj) are 

R P \\rvK a p]p^ + R P p\ TV ^a(i\\p = . (34) 

The extra algebraic relations ()32|) can be treated in the manner explained in section 17.21 
by defining the tensor 

and imposing P a pxp = as initial conditions. This can be shown to be preserved by the 
evolution, see [TJ |§] . 

Similarly, the Einstein equations themselves can be treated as initial conditions by the 
method of subsection 17.31 whenever the explicit form of the energy-momentum tensor in 
terms of the gravitating physical fields is given. (This has been done in 4 dimensions for 
some cases in |TT1 ITS].) I am going to illustrate the method by considering the case in 
which there is only a massless scalar field <fi on the spacetime. 

The energy-momentum tensor of the scalar field is given by ()26|) . hence the needed 
choice for jpxp is, using the notation introduced at the end of subsection 17.51 

This suggests that the proper thing to do is going one step higher for the scalar field 
defining 

4>tiu = , 4>tiu = 4>{tiu) 

whose field equations are easily determined using the method of subsection 17.51 

V[a<^]„ = ~4> P R p uXn , V p </v = -4> p R P u ■ 

The Einstein field equations read simply 

8nG 1 1 2 a 

Rpu ~[-<Pn<Pu + ——^^9^ = . 

Therefore, the system to be considered is a mixed system, for the unknowns (K a/ 3Xp, 4>p,v, 4>p) 
where K a/3X p is a double (2,2)-form, is a double symmetric (l,l)-form and M is a one- 
form. The equations are (on a given background) 

1 , „„ „„ , , ( 2 . 87rG 

pv = <Pl 

16itG 



v M <^ = (j)pu , v [A </v]„ = -^ P R p v\p , v p <frpv = <\>v \—^7^ — ~^^ P<i 



VfrKaftXp = , VKppXp, 



c 4 
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This must be supplemented with the extra algebraic relations (|3*4*)l together with 



4> p p = 0, K[c#\\n = 0, K \^v)p + -^4>p4>u + ^7^ A ^ = 

all of which can be considered as initial conditions for the data on an initial spacelike 
hypersurface E. Of course, this data must be restricted to satisfy the constraint equations 
provided by (|1()H11|) : 

N [x (V^ - (f)^) = 0, N [T V x ^]u + \<P P R p u[XpN T] = 0, 
iV[rV 7 K a/3]v = 0, {v p K pf3X , - ±^0^) = 

where N is the normal to S. This system of equations is directly hyperbolizable as 
shown in previous sections by using superenergy tensors. Here, apart from the Bel tensor 
(J30J) mentioned above, one has to use the superenergy tensor of the scalar field (jHSj and 
references therein): 

T a p\n {0[i][i]} = 20 Q ( A At ) /3 - g a ^x p (p^p - g\p(f) a p <frpp + -g a pgx^ p <t>cr P ■ 

Observe that both needed superenergy tensors have the same number of indices and the 
same symmetry properties. They can in fact be added into one single, mixed, superenergy 
tensor 

= T af3\p {K[2] [2]} + T a p Xfl {</» [i][i]} 

whose total superenergy density relative to a timelike unit vector u, and the total w(t) 
of section El is the sum of the separate corresponding quantities for each of the two 
fields involved. It is very interesting to remark that the inequalities computed in section 
El can then be easily estimated by using the divergence of the total tensor above [HE] . 
Furthermore, in cases with Killing vectors, the corresponding total, mixed, superenergy 
currents are divergence free, and provide conserved quantities, see subsection 7.2 in [53] , 
A similar procedure can be used to couple the Bianchi equations to other matter fields, 
such as the electromagnetic one. In that case, one must use the symmetrized superenergy 
tensor associated to the covariant derivative of the Maxwell 2-form, which is the Chevreton 
tensor ^J|S3]. As a matter of fact, one can go to higher-order derivatives of all the fields 
involved (the Riemann tensor, the scalar field, the electromagentic 2-form, etcetera) and 
find hyperbolizations of the new variables by using the different (super) ^-energy tensors, 
for all natural numbers k, introduced in [35] . 
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